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Summary. — In this talk I will review the "state of the art" of the calculations of 
the nuclear matrix elements (NME) of the neutrinoless double beta decays (0u(3/3) 
for the nuclei 48 Ca, 76 Ge, 82 Se, 124 Sn, 128 Te, 130 Te and 136 Xe in the framework of 
the Interacting Shell Model (ISM), and compare them with the NME's obtained 
using the Quasi-particle RPA approach (QRPA). I will also discuss the effect of the 
competition between the pairing and quadrupole correlations in the value of these 
NME's. In particular I will show that, as the difference in deformation between 
parent and grand daughter grows, the NME's of both the neutrinoless and the two 
neutrino modes decrease rapidly. 



1. Introduction 

The discovery of neutrino oscillations in recent experiments at Super-Kamiokande 
PQ, SNO [5] and KamLAND [3] has changed the old conception of neutrinos by proving 
that they are massive particles. According to the origin of their mass, neutrinos can be 
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either Dirac or Majorana particles, the latter case being particularly interesting since it 
would imply an extension to the standard model of electroweak interactions, and, being 
neutrinos their own antiparticles in this scenario, lepton number conservation would be 
broken. Besides, it happens that the best way to detect one of these violating processes 
and consequently establish the Majorana character of the neutrinos would be detection 
of the neutrinoless double beta decay [Quf}f3). 

Double beta decay is a very slow weak process. It takes place between two even-even 
isobars when the single beta decay is energetically forbidden or hindered by large spin 
difference. Two neutrinos beta decay is a second order weak process — the reason of its 
low rate — , and has been measured in a few nuclei. The §v(3(3 decay is analog but needs 
neutrinos to be Majorana particles. With the exception of one unconfirmed claim [H[5], 
it has never been observed, and currently there is a number of experiments either taking 
place [BJ El [5] or expected for the near future — see e.g. ref. [9] — devoted to detect this 
processes and to set up firmly the nature of neutrinos. 

Furthermore, 0v/3{3 decay is also sensitive to the absolute scale of neutrino mass, and 
hence to the mass hierarchy — at present, only the difference between different mass 
eigenstates is known. Since the half-life of the decay is determined, together with the 
masses, by the nuclear matrix element (NME) for this process, the knowledge of these 
NME's is essential to predict the most favorable decays and, once detection is achieved, 
to settle the neutrino mass scale and hierarchy. 

Two different and complementary methods arc mainly used to calculate NME's for 
Oiy/3/3 decays. One is the family of the quasiparticle random-phase approximation (QRPA). 
This method has been used by different groups and a variety of techniques is employed, 
with results for most of the possible emitters [10l [Til H2] ■ This work concerns to the 
alternative, the interacting shell model (ISM) [15] , 

In previous works [TUfTS], the NME's for the 0^/3/3 decay were calculated taking into 
account only the dominant terms of the nucleon current. However, in ref. 16J it was 
noted that the higher order contributions to the current (HOC) are not negligible and it 
was claimed that they could reduce up to 20%-30% the final NME's. Subsequently, other 
QRPA calculations jTTJ [TS] have also taken into account these terms, although resulting 
in a somewhat smaller correction. These additional nucleon current contributions have 
been recently included for the first time in the ISM framework [TS]. In addition, the 
short range correlations (SRC), are now modeled cither by the Jastrow prescription or 
by the UCOM method pb]. 

2. ISM vs QRPA Nuclear Matrix Elements 

The expression for the half- life of the 0v/3j3 decay can be written as PU H2 : 



where (m„) = U^rrik is the averaged neutrino mass, a combination of the neutrino 
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masses m,k due to the neutrino mixing matrix U — as we see, the neutrino mass scale is 
directly related to the decay rate — and Gqi is a kinematic factor — dependent on the 
charge, mass and available energy of the process. M ""* is the NME object of study in 
this work. 

The kinematic factor Goi depends on the value of the coupling constant qa- Therefore 
we have to take this into account when comparing the values of NME's obtained with 
different qa values. In these cases we will use a NME modified as: 

(2) M' 0v W = (jt^) 2 M 0l,f313 

These M' 0l//3/3 's are directly comparable between them no matter which was the 
value of qa employed in their calculation, since they share a common Gqi factor — that 
°f 9a = 1-25. In this sense, the translation of M /0l/ ^'s into half-lives is transparent. 
The QRPA results obtained with different gA values are already expressed in this way 
by the authors of refs. [UJ HI] while the results of refs. (TTJ, HHl [H] will be translated by 
us into the above form when compared with other results. 

We have calculated the ISM NME's both taking the UCOM and Jastrow ansatzs for 
the short range correlations. The former use the correlator of the ST = 01 channel [25] . 
throughout the calculation. The correlator of the other important — even — channel is 
very similar to this one, and it should not make much difference on this result. 

In figure[T]the ISM and QRPA results for the NME's are compared within this UCOM 
treatment of the SRC. The same figure but considering Jastrow type SRC was shown 
in ref. [19], but, inadvertently, the QRPA results from refs. [T71 HH1 121] obtained with 
9,4=1.0 where not transformed properly. This is corrected in figure[2j By comparing both 
figures, it is confirmed that there is a common trend; when the nuclei that participate in 
the decay have a low level of quadrupole correlations, as in the decays of 124 Sn and 136 Xe, 
both approaches agree. The QRPA in a spherical basis seems not to be able to capture 
the totality of the quadrupole correlations when they are strong. As these correlations 
tend to reduce the NME's, the QRPA produces NME's that are too large in 76 Ge, 82 Se, 
128 Te, and 130 Te. For both ISM and QRPA the only net effect of UCOM is an increase 
of the Jastrow results of about 20%. 

3. The Influence of Deformation in the NME's 

An important issue regarding Ou/3/3 decay is the role of pairing and deformation. It 
has recently been discussed in ref. [TpJ that the pairing interaction favors the Qvf3j3 decay 
and that, consequently, truncations in seniority, not including the anti-pairing- like effect 
of the missing uncoupled pairs, tend to overestimate the value of the NME's. On the 
other hand, the NME is also reduced when the parent and grand-daughter nuclei have 
different deformations |26] . Thus, we have studied the interplay between both pairing and 
deformation, this is, to which extent a wave function in the laboratory frame, truncated 
in seniority, can capture the correlations induced by the quadrupole-quadrupole part of 
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Fig. 1. - The neutrinoless double beta decay M' °"' 3 ' 3 's for ISM and QRPA calculations treating 
the SRC with the UCOM approach. Tu07 QRPA results from ref. [23] and Jy07 results from 
refs. [TanS. 



the nuclear interaction, and its eventual influence in Ou/3/3 NME's. There is an extra 
motivation to pursue this study; the possibility to carry on an experiment with 150 Nd, 
which is a well deformed nuclei, decaying into 150 Sm which is a much less deformed one. 

To study the interplay between pairing, seniority truncations, and quadrupole corre- 
lations we need first to decide how to measure the quadrupole correlations of the ground 
state. Our choice is to refer to non energy-weighted sum rule: 



(3) <Q 2 ) = £ |<2+| Q|0+)| 2 

i 

The operator Q represents the mass quadrupole. No effective "nuclear" charges are 
included. Using 82 Kr as our test bench, we proceed to compute (Q 2 ), first with our 
standing effective interaction and different seniority truncations (s m means the maximum 
seniority allowed in the wave functions of parent and grand daughter nuclei) . The results 
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Fig. 2. - Same as fig[T]but with Jastrow type SRC. Tu07 QRPA results from ref. [12] and Jy08 
results from ref. [24] , 



are drawn in Figure [3] as the black circles labeled A = 0. We can see that at s m = 4 
— roughly, the implicit level of seniority truncation in the spherical QRPA — some 70% 
of the full quadrupole correlations are incorporated in the wave function. We would like 
to know how this behavior evolves when more correlations are enforced in the system. 
For this we recalculate the ground state of 82 Kr with a new hamiltonian that consists of 
the standing one plus a quadrupole-quadrupole term XQ ■ Q, whose effect will be gauged 
by its influence in the sum rule. To have an idea of the relevant range of values of (Q 2 ) in 
this nucleus and valence space, we have gone to the limit of pure quadrupole-quadrupole 
interaction with degenerate single particle energies, getting (Q 2 ) ~ 4500 fm 4 . The results 
for A = 1 and A = 2 are also shown in Figure [3] (A = 1 corresponds to X qq — 0.025 in 
Figs. 6-9 and \ qq = 1 in Fig. 10 to X qq =0.1. It is evident in the figure that, as we try 
to increase the correlations, the s m — 4 truncation becomes more and more ineffective. 
For A = 1, only 57% of the full correlations are present, and for A = 2 only 50%. The 
situation is different for 82 Se: while for A = the values of (Q 2 ) as a function of seniority 
are similar, albeit a bit smaller than the 82 Kr ones, for A = 1 and A = 2 there is scarcely 
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Fig. 3. - Quadrupole correlations in the ground state of 82 Kr as a function of the amount 
of quadrupole-quadrupole interaction XQ ■ Q added to the hamiltonian and of the maximum 
seniority s m permitted in the wave functions 



any increase of the ground state correlations. This means also that, as we increase A, 
the "deformation" of 82 Kr grows, whereas that of 82 Se remains constant. 

This behavior offers us the opportunity of exploring the effect of the difference in 
deformation between parent and grand daughter in the 0v(3(3 NME's. To this goal, 
we have computed the Gamow- Teller matrix element for different values of A — the 
amount of extra quadrupole-quadrupole interaction — and s m — the maximum seniority 
allowed in the wave function — . The results are gathered in Figure 2J For s m = 0, we 
observe that the Gamow- Teller matrix element grows as a function of A. This may seem 
paradoxical, but is not, because at this seniority truncation, the only effect of adding 
more quadrupole-quadrupole interaction is to augment the pairing content of the wave 
functions, thus increasing M GT . At s m = 4, M GT remains constant as a function of A, 
meaning that the minor increase of the correlations of 82 Kr, that we have shown in Figure 
is barely enough to compensate the increase of M GT at s m — 0. On the contrary, the 
full space results are sensitive to the difference in deformation — or, to be more precise, 
to the difference in the level of quadrupole correlations in the ground state — between 
parent and grand daughter. The effect goes in the direction of reducing the value of 
M GT . In the A = 82 decay, doubling the quadrupole correlations in 82 Kr, roughly halves 
M GT . 
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Fig. 4. — Se — v Kr Gamow- Teller matrix element, M , as a function of the maximum 
seniority of the wave functions, for different values of the strength of the extra quadrupole- 
quadrupole interaction. 



We can go much further in the exploration of the deformation effects, calculating the 
NME of the decay for initial and final states computed with different amounts of sup- 
plementary quadrupole-quadrupole interaction. As before , we measure the quadrupolc 
correlations by means of the mass quadrupole sum rule. The results of this search are 
plotted in Fig. [Sj We observe that the NME decreases almost linearly as the difference 
of the sum rules for the final and initial states increases. In fact, the maximum values 
of the NME are reached when this difference is close to zero. On the contrary for large 
differences the NME can be extremely quenched. Notice that the values in the figure 
range between 0.07 and 2.7, a factor of 40 span. The A=0 value is 2.18, indicating that 
the difference in quadrupole correlations between 82 Kr and 82 Se is not very large. 

4. Of (Unphysical) Mirror Decays: A Case Study 

We have also studied the transitions between mirror nuclei in order to have a clearer 
view of the role of deformations in the NME's. These transitions have the peculiarity 
that the wave functions of the initial and final nuclei are identical (provided Coulomb 
effects arc neglected) and consequently the interplay of the Qv(3j3 operator and of the 
nuclear wave functions in the NME may be easier to understand. 

We have studied four parent nuclei with six valence protons and four valence neutrons, 
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Fig. 5. - 82 Se — » 82 Kr NME, M 0l/ , as a function of the difference between the mass quadrupole 
sum rule between 82 Kr and 82 Se for a large number of different values of the strength of the 
added quadrupole quadrupole interaction. 



26 Mg, 50 Cr, 66 Ge and 110 Xe, decaying into the four grand daughter nuclei with four 
valence protons and six valence neutrons, 26 Si, 50 Fe, 66 Se and 110 Ba. The valence spaces 
considered are the sd-shell, the p/-shell, r 3 g (IP3/2, lpi/2> 0f 5 / 2 , 0g g / 2 ) and r4h (0g 7 / 2 , 
ld 5/2 , ld 3/2 , 2si /2 , 0h u/2 ). The starting interactions are USD [57], KB3 [23], GCN28.50 
and GCN50.82 [29] . The deformation of the nuclei is modified as in the previous section, 
and it is quantified also in the same way. 

The results for A = 66 in the case of equally deformed initial and final nuclei, are 
shown on Fig. [6] There we see that, as the nuclei become more deformed, the NME 
and the pairing content of the wave function get smaller, while the quadrupole sum rule 
grows. All these changes are nearly linear for reasonable deformations and then the 
saturation is approached more smoothly. Note that the purely quadrupole interaction 
(X qq ~>-oo limit) gives a NME which is about a half of the value obtaind with no additional 
quadrupole. 

Fig. [7J shows the same quantities than Fig. [5] but now only the final nucleus has been 
artificially deformed by adding an extra quadrupole-quadupole term. In addition, the 
overlap between initial and final wave functions has been included. We see that now, the 
reduction of the NME is more pronounced and, what is more interesting, that it follows 
closely the overlap between wave functions. This means that, if we write the final wave 
function as: \^) = a |^o) + b\^ qq ), the 0vf3f3 operator only connects the parts of the 
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Fig. 6. 



3 Ge 



3 Se NME, M ", as a function of the strength of the added quadrupole 



quadrupole interaction. Equally deformed case; the same amount of extra QQ interaction is 
added to 66 Ge and 66 Se. On the right hand y axis the pairing and quadrupole sum rules are 
represented, normalized so that their maximum value is 1. Note the change of scale in the x 
axis at A = 0.2. 



wave functions that have the same deformation among themselves. 

The behavior of the NME's with respect to the difference of deformation between 
parent and grand daughter is common to all the other transitions between mirror nuclei 
that we have studied. Therefore we can submit that this is a robust result. However, 
when we consider the transitions between equally deformed nuclei, the evolution of the 
NME's with the deformation that we have found in A=66 is only shared by the A=110 
case. When the valence space is a full major oscillator shell — A = 26 and A = 50 — 
the situation is quite different. Indeed, what is observed is that the NME does not 
decrease for moderate values of A but remains rather constant until a point — with large 
deformation — where its value increases significantly — up to 50% — . This is due to the 
fact that, at this point, the major contribution to the NME ceases to come only form the 
decay of pairs coupled to J = 0, since other values like J = 2,4,6, which usually have 
a contribution to the NME contrary to that of J = 0, reverse sign and grow until being 
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comparable with this contribution, thus resulting in this notorious rise of the NME. In 
the sd-shell and p/-shell cases, the A gg ^oo limit is equivalent to Elliott's SU(3) limit, 
and, the fact that both the initial and final nuclei should belong to the same irrep of 
SU(3) may be the reason of the increase of the NME, but, for the moment we have not 
found a formal explanation. 

5. 2v (Unphysical) Mirror Decays: A Case Study 

Very similar conclusions may be reached for the effect of deformation in 2vf3f3 decay. 
For instance, the mirror nuclei diagonal and non-diagonal NME's are represented as in 
the 0vPf3 case in Figs. [5] and [9] for the A = 66 transition. We see that the figures 
resemble very much that of the previous section, with the only exception that, in the 
equally deformed, case, the lowering of the NME due to deformation is more pronounced. 

If we look to a non mirror — but again fictitious — transition, for instance that of 
48 Ti —* 48 Cr, we get the results shown in Fig. [101 If we compare with those of ref. [26] . 



Deformation and the Nuclear Matrix Elements of the Neutrinoless /3/3 Decay 



11 




which are again the equivalent ones for the Of/3/3 transition, we see that there are not 
substantial changes. In this sense, deformation seems to affect similarly to Qv(3fj and 
2v(i(3 decays. 



6. Summary 

After a brief discussion of the "state of the art" results for the nuclear matrix elements 
of the neutrinoless double beta decay in the context of the Interacting Shell Model and 
of the Quasiparticle Random Phase Approximation, we have analyzed the role of the 
pairing correlations and the deformation in the NME's, concluding that seniority trun- 
cations are less reliable when the quadrupole correlations are large. Since the NME's are 
reduced when the deformation of parent and grand daughter is different, a bad treat- 
ment of the quadrupole correlations can lead to an artificial enhancement of the NME's 
in the transitions among nuclei with unequal deformations, or in the cases of nuclei with 
different -and large- amounts of quadrupole correlations. 
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Fig. 9. - The same as the previous figure, but now the only nuclei calculated with additional 
quadrupole interaction is the final one. The normalized overlap between initial and final states 
is also included. 
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